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1. Introduction 

Enumeration of Hamiltonian walks (HWs), i.e. self-avoiding walks (SAWs) that visit 
every site of a given lattice, is a classic problem in graph theory, but it also has an 
important role in the study of the configurational statistics of polymers. HWs are 
used to model collapsed polymers [JJ, polymer melting 0111], as well as protein folding 
jU |3] . The number of all possible HWs on a lattice is related to the configurational 
entropy of a collapsed polymer system, and also to the optimal solutions to the 
traveling salesman problem jHJ Ej . Enumeration of HWs, closed or open, is a difficult 
combinatorial problem, which has been exactly solved only for few lattices, namely, the 
two-dimensional Manhattan oriented square lattice jHUH], the two-dimensional ice lattice 
[TU] . the two-dimensional hexagonal lattice [IHH2], the Sierpinski gasket fractal jHj, and 
the 4-simplex fractal ^3]- The number of HWs has also been calculated numerically 
for various lattices by means of direct enumeration [TH ITS] , transfer matrix methods 
[TBHmHH], and Monte Carlo estimates [IHlEIIj. The field theory representation for this 
problem was introduced in [3TJ, and further developed in [22] and [2l!j . 

The purpose of this paper is to understand the topological properties of Hamiltonian 
walks on several families of fractal lattices. In order to do this we study the asymptotic 
behavior of the number of closed HWs for a large number of vertices N. This analysis 
yields the values of the so-called connectivity constant u which has the physical meaning 
of the average number of steps available to the walker having already completed a large 
number of steps. It also provides insight into the spatial distribution of walks present 
at large N; this can be related to detailed studies of knot delocalisation in Refs. [2I| . 

The number is for homogeneous lattices with N 3> 1 expected to take the form 

C N ~u N ns m N a . (1.1) 

Here uo is the connectivity constant and the term with /is represents a surface correction 
(fis < 1), with a = (d — l)/d (d being the dimensionality of the lattice). This differs 
from the ordinary SAW case, where no surface term fJ,s N " is expected, i.e. the number 
of SAWs of length iV behaves as fi N N a for large N. Furthermore, the exponent a 
is universal in the SAW case, i.e. it depends only on the dimensionality of the lattice, 
whereas for HWs it may depend on other, not yet identified characteristics of the lattice. 
To the lowest approximation in Eq. (jl.lj) Cn ~ w , and the connectivity constant can 
be defined as 

l nu= lim i^k. (1.2) 

N^oo N 

For a better understanding of these problems it is helpful to study HWs on fractal 
lattices. As was first recognized by Bradley ^3], the self-similarity of fractal lattices is a 
useful tool for the exact and computationally fast iterative generation and enumeration 
of all HWs on an unlimitedly large corresponding fractal structure. In this paper we 
extend Bradley's algorithm to two- and three-dimensional Sierpinski fractal families, 
as well as n-simplex fractals with n > 4. We consider several families of fractals in 
order to compare the obtained results and be able to draw more general conclusions 
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about the character of the Hamiltonian walks on different classes of lattices. In the 
case of the two-dimensional Sierpinski fractal family an exact closed form result for the 
connectivity constant is obtained due to the simple form of the recursion relations for 
the numbers of HWs. For the three-dimensional Sierpinski fractals which can model 
physically more frequently encountered systems a numerical approach is necessary. The 
study of asymptotics of the number of HWs shows that the surface term in Eq. (jl.ljl 
only appears for fractals on which the collapse transition from the polymer coil to the 
globule phase is possible. This is the same class of lattices for which delocalised HWs 
dominate over localised ones for large N. 

The paper is organized as follows. In section 2 we describe the two-dimensional 
Sierpinski fractal family and obtain exact recursion relations for the number of HWs, as 
well as the closed exact formula for to. Recursion relations for HWs on three-dimensional 
Sierpinski fractals are given and analyzed in section 3. A similar method for analyzing 
HWs on 5- and 6-simplex lattices is presented in section 4. Finally, in section 5 we 
discuss all our findings and related results obtained by other authors. 



2. Hamiltonian walks on two— dimensional Sierpinski fractals 

We begin by defining the two-dimensional (2d) Sierpinski fractal family. Each member 
of the 2d SF family (labeled by 6) can be constructed recursively, starting with an 
equilateral triangle that contains 6 2 smaller equilateral triangles (generator G\(b)). The 
subsequent fractal stages are constructed self-similarly, by replacing each of the 6(6+1) /2 
upward-oriented small triangles of the initial generator by a new generator. To obtain 
the /th-stage fractal lattice Gf(b), which we shall call the /th order generator, this 
process of construction has to be repeated I — 1 times, and the complete fractal is 
obtained in the limit I —>■ oo. It is easy to see that, for any 2d SF, each closed HW on 
the (I + l)th generator is comprised of HWs which enter and exit Zth order generators. 
Let C\ be the number of closed HWs on the /th order generator, whereas h\ and g\ are 
the numbers of HWs which enter the /th order generator at one vertex, and leave it at 
the other, with or without visiting its third vertex, respectively. Then, it can be shown 
that a simple relation 

C l+l = Bh«g? (2.1) 

is valid for / > 1 (see Appendix A). Here B is a constant that depends only on SF 
parameter b, whereas exponents a and (3 are equal to 

a = 6+ 1, „=£±^. (2.2) 

For instance, the explicit form of the relation (|2.1j) for b = 2 is Ci + \ = hf and for 6 = 3: 
Ci+i = 3hfgf, which is illustrated in Figure [TJ One can also show (Appendix A) that 
numbers hi and gi of open HWs obey the following closed set of recursion relations 

h l+1 = Ahfgf, g l+1 = ANT 1 it 1 , ( 2 -3) 
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(a) 






(b) 



Figure 1. The (/ + l)th order generators G;+i(2) (a) and Gj+i(3) (b) for 6 = 2 and 
b = 3 two-dimensional Sierpinski fractals, with the possible closed Hamiltonian walks 
configurations depicted. The small up-oriented triangles are the respective Ith order 
generators, and the lines that traverse them represent the open Hamiltonian walks. 











Figure 2. Schematic representation of the recursion relation for the number of 
open Hamiltonian walks (2.3) of h— type on b = 3 two-dimensional Sierpinski fractal 
structures. 



for all / > 1, where A is again a constant, different for every 6, and 

6(6-1) 



x 



y 



(2.4) 



For 6 = 2 these relations have the form: hi + i = 2hfgi, gi + i = 2higf and for 6 = 3: 
hi+\ = 8hfgf, gi + i = 8hfgf (see Figure |2J). From the relations (|2.3|) it follows 
straightforwardly that 



H = iL = k 
hi h x 



(2.5) 



for any I > 1, so that from J2DJ and (Q one gets C l+1 = BK p hf >+1)/2 . Since the 
number N\ of sites on the Zth order generator satisfies the recursion relation 

6(6 + 1) 



N, 



-M " (b 2 ~ 1), 



(2.6) 



according to f)1.2|) it then follows that 



\iiuj = lim 



InC, 



l+i 



lim 



In hi 



l^co Ni + \ i^oo Ni 



Hamiltonian walks on fractals 



5 



Table 1. Values of the number A of Hamiltonian configurations relevant for the 
recursion relations (2.3), numbers hi and g\ of open HWs on the generators of two- 
dimensional Sierpinski fractals, connectivity constant oj for HWs, and connectivity 
constant /i for SAWs (obtained via RG method - /i RG and numerically estimated n nnm 



in 


28]), for 2 < b < I 










b 


A 


hi 




u> 


|U RG (/i num ) 


2 


2 


2 


3 


1.31798 


2.288(2.282 ±0.007) 


3 


8 


10 


11 


1.39157 


2.491(2.49 ±0.02) 


4 


40 


92 


112 


1.46186 


2.656(2.686 ±0.004) 


5 


360 


1 852 


2 286 


1.52155 


2.791(2.82 ±0.01) 


G 


3 872 


78 032 


94 696 


1.56895 


2.904 (2.92 ±0.02) 


7 


62 848 


6 846 876 


8 320 626 


1.61011 


3.005(2.99 ±0.05) 


8 


1 287 840 1 255 156 712 


1 527 633 172 


1.64528 


— (3.13 ±0.07) 



On the other hand, from (Q, (JUj), and (Q one has h l+1 = AKyh h[b+1)/2 , i.e. \nh t 
satisfies the difference equation lnh[ + i = lnAK y + b( - b ^ In hi, whose solution is 

+ 1 — [6(6 + 1)/2] v ; V 2 / 
From this equation, together with the explicit expression for the number of sites 

6 + 4 /K6 + 1)V +1 + £M 
+ 6 + 2 V 2 y 6 + 2' v ' 

which follows directly from the relation (|2.fi|) . one can derive the general form 

4 4 2 

w = A b{b + 4)(6= - 1) h Kb + 1)(6 + 4) ^(6 + 1)(6 + 4) (2 Q) 

of the connectivity constant (jl.2j) for HWs on two-dimensional SFs. Consequently, in 
order to calculate u for any particular 2d SF, one should find the numbers h\ and gi of 
open HWs on the generator, and the number A of all Hamiltonian configurations which 
are relevant for recursion relations ()2.3J) . In Table 1 we present these numbers, together 
with the values of lu, for 2 < 6 < 8. As one can see, for 6 = 2, and 6 = 3 the numbers 
gi, hi, and A are small and can be directly enumerated, whereas for larger values of 6 
they quickly increase, so that enumeration should be computerized (calculation of the 
numbers A, gi, and hi required 13 minutes for 6 = 7 case, and about 100 hours for 6 = 8, 
both on a computer with a processor MIPS R10000, Rev 2.6 on 180 MHz). One should 
mention that the connectivity constant for the Sierpinski gasket (6 = 2) has already 
been calculated in a different way by Bradley 

Combining equations ()2.7j) . ()2.8j) and ()2.9j) it is not difficult to see that hi = Gu Nl , 
where G depends only on 6, and correspondingly C\ ~ to l . Comparing with ()1.1|) one 
can conclude that neither surface nor power correction terms are present in the scaling 
form for the number of closed HWs on two-dimensional Sierpinski fractals. 
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Comparison with the self- avoiding walk case 

It is interesting to compare the value of the connectivity constant obtained for the case 
of Hamiltonian walks to that corresponding to all possible self-avoiding configurations 
on 2d SFs. Our algorithm for enumerating HWs is easily adjusted for that purpose. 
By means of exact renormalization group (RG) approach |2"5| EE] , these configurations 
can be used for calculating the connectivity constant \x for ordinary SAWs, which was 
done earlier only for the 6 = 2 case [2E] . Here we extend an exact RG calculation of the 
connectivity constant /x for any 6. 

The connectivity constant \i for the SAW model is equal to \x = 
limAr_ >00 (cAr + i/cAr) = liniAr^oo (pn+i/pn) , where cn (pn) is the average number of dis- 
tinct open (closed) n-step SAWs. In order to calculate /i within the exact RG approach, 
one should introduce two generating functions and B[ : 

B« = B? = J2B® N x N , 

N N 

where x is the statistical weight of each step of the SAW (fugacity), whereas B$ (Bi n) 
is the number of SAWs which enter the Ith order generator Gf(b) at one vertex, and 
leave it at the second, without (with) visiting the third one. For every 2d SF lattice 
functions B and Bi obey recursion relations of the following form 

b {i+i) = e {B^y (B?y , B r i) = $>„■(&) o«y (B?y , ( 2 .io) 

where fij(b) and gij(b) are coefficients that do not depend on I, but do depend on the 
fractal parameter b. The initial conditions are B^ = x, B^ = x 2 and the connectivity 
constant /i is equal to 1/x*, where x* is the value of the fugacity for which one approaches 
the fixed point (B*,B*) of (|2.10p . after a large (infinite) number of iterations. 
The explicit RG recursion relations for 2d SFs with 6 = 2 are 

B' = B 2 + B 3 + 2BB X + 2B 2 Bi + B{, B[ = B 2 B l + 2BB 2 

and for 6 = 3: 

B' = B 3 + 3B 4 + B 5 + 2B 6 + 3B 2 B 1 + 12B 3 B 1 + AB A B 1 + %B b B l + 3BB 2 + lQB 2 Bj 

+ 5B 3 B 2 + 8B 4 B 2 + +B\ + 8BBf + 2B 2 Bf + Bf 
B[ = B A B l + 2B 5 B 1 + AB 3 B\ + 8B 4 B 2 + 5B 2 Bf + 8B 3 B 3 + 2BB\ 

whereas relations for 6 = 4 and 5 are given in Appendix B. For 6 = 6 and 7 RG relations 
are too cumbersome to be quoted here, but they are available upon request. Underlined 
terms in quoted RG relations correspond to the Hamiltonian configurations, as one can 
check by direct comparison with relations (|2.Hj) : replacing gi and hi in (|2.3jl by B and 
£>i, respectively, one obtains underlined terms in the RG relations (values of coefficient 
A are given in table 1). It is obvious that for larger 6 values the number of SAW 
configurations is much larger than the number of Hamiltonian configurations, and for 
that reason we were not able to enumerate all SAW configurations on 2d SF beyond 
6 = 7. 
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For all b considered here recursive relations (|2.1Uj) have only one nontrivial fixed 
point (B%,0), where B\ lies in the interval < B\ < 1. One can check that by setting 
B\ = in (j2.1()j) RG equations used in |2Zj for calculating critical exponent v (connected 
with the mean end-to-end distance) for SAWs on 2d SFs are recovered. Of course, SAW 
model treated in |2Zj is slightly different (each unit triangle within the fractal can be 
traversed only along one side) from the usual one, treated here, but both of them belong 
to the same universality class, i.e. the critical exponent v is equal for both considered 
SAW models. This is not the case with the connectivity constant, which is a nonuniversal 
quantity, so an extra RG parameter B 1 had to be introduced. 

The final RG results [i for the SAW connectivity constant are given in the last 
column of Table 1, together with the corresponding values ^ num obtained in [2E| using 
a graph counting technique. As one can expect, the connectivity constant /i for SAW 
is larger than uj for HW model, for every considered SF, since the physical meaning of 
the connectivity constant is the average number of steps available to the walker after 
N steps completed, for large N. One can also see that the values of /i obtained by two 
different methods for 6 = 4 and 5 are not in good agreement. The RG method applied 
here is exact, implying that numerical estimations within the graph counting technique 
used in j2H] were not accurate enough. 

3. Hamiltonian walks on three— dimensional Sierpinski fractals 

We proceed by analyzing the properties of HWs on the three-dimensional (3d) SF 
family. A 3d SF can be constructed recursively as in the 2d case, the only difference 
being that the generator G\{b) of the fractal with the parameter b is no longer a triangle, 
but a tetrahedron that contains b(b+ l)(6 + 2)/6 upward oriented smaller tetrahedrons. 
Consequently, one should observe four types of open HWs that traverse the Z-th order 
generator in order to obtain the overall number of closed HWs on 3d SF. The first 
three types are HWs which enter the generator at one vertex and leave it at another, 
meanwhile 

• visiting the third, but not the fourth vertex (g-type), 

• visiting both the third and the fourth vertex (/i-type), 

• visiting neither the third nor the fourth vertex (z-type). 

HWs of the fourth possible type (j-type) consist of two self-avoiding branches, and 
correspond to the walks that enter the generator and leave it without visiting the 
remaining two vertices, then, after visiting other parts of the lattice, enter the same 
generator again at the third corner vertex and finally leave it at the fourth corner 
vertex. Examples of these types of walks are sketched in figure 3. In principle, it is 
possible to establish a closed set of recursion relations for the numbers gi,hi,ii, and ji 
of the corresponding walks for any 3d SF in the following form 



J2 g(n,m,k)g?hrifji 



6(ft+l)(ft+2) 



(n+m+fc) 



9i+i = 
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i>Q+l)Q+2) , , , 

/i m = £ W(n, m, fc)lWfo 6 > 

n,m,fc 
n,m,k 

b(b+l)(b+2) , , s 

j m = £ jT(n,m,%r^ 6 ( \ (3-1) 

where Q(n,m,k), TC(n,m,k), T(n,m,k), and J(n,m,k) are the numbers of open 
Hamiltonian configurations of the corresponding types, with n branches of the g— type, 
m branches of the fo-type, branches of the i-type, and [b(b + 1) {b + 2)/6 — (n + m + k)} 
branches of the j-type. The number Ci + i of all closed HWs within the (/ + l)th order 
generator for any b is equal to 

b(b+l)(&+2) , 

C w = Yl B{n,m,k)g^hTitji e ( + + ) , 

n,m,k 

where B(n, m, k) is the number of all closed Hamiltonian configurations with n branches 
of the g-type, m branches of the h— type, k branches of the z— type, and [b(b + l)(b + 
2)/6 — (n + m + k)] branches of the j-type. As one can see, the recursion relations (j3.1|) 
for the number of open HWs on 3d Sierpinski fractals are much more complicated than 
the corresponding equations ()2.3|) for 2d SFs. Consequently, it is not possible to find 
explicit expression, similar to ()2.9|) . for the connectivity constant uj. Instead, one should 
perform a numerical analysis of the recursion relations (J3.1j) in order to find the value 
of uj. We shall demonstrate the method on the particular case of the 6 = 2 fractal. 

By computer enumeration of the possible HW configurations within the (I + l)th- 
order generator Gf +1 (2) of the b = 2 3d Sierpinski fractal, we found the following 
recursion relations: 

g ' = Q g 2 f + Ag 3 j + 2g A + Ylifh + 24igjh + 2Aig 2 h + 8i 2 h 2 , (3.2) 

h! = 24fhg + 16h 2 ij + 16hg 3 + 32h 2 gi + 24g 2 hj , (3.3) 

i' = \2igj 2 + 12^ 2 j + 8ig 3 + 8i 2 jh + 16i 2 gh , (3.4) 

f = 8i 2 h 2 + A8igjh + 22j 4 + 2g 4 + 8g 3 j + 2Aig 2 h , (3.5) 

where we have used the prime symbol as a superscript for the numbers of HWs on the 
(/ + l)th-order generator Gf +l (2) and no indices for the Zth-order numbers. From the 
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definition (|1.2j) of the connectivity constant and the formula for the number C/ + i of 
closed HWs within Gf +1 (2): 

C l+1 = 16gfh? , (3.6) 

it then follows that 

lnC m 1 ln# 1 ln/i, 

mo; = hm — — = - lim — — — — hm — — , (3.7) 

where Ni = 2(4' + 1) is the number of sites in Gf(2). On the other hand, from the 
recursion relation ()3.2|) for the numbers gi one obtains 

Inflfj+i 1 ln# 1 lnj; 
hm = - hm 1 — iim 

i^oo Ni + i 2 i^oo Ni 2 i^oo Ni 



+ - hm -r-rln 6 + Axi + 2xj + 2Ayi zi -\ = 1 1 ^— (3.8) 

2 i^oo 4' +i \ xf xi xf J 

where x\ = gi/ji, yi = hi/ji and z\ = ii/ji- In a similar way, from ()3.3)1 it follows that: 

ln/ij+i 1 ln/i; 1 ln^j 1 lnj; 
hm = - lim 1 — hm 1 — lim 

i^oo Ni + i 4 i^oo Ni 4 i^oo Ni 2 i^oo Ni 



1 , : _ 1 ,.. o fn , n .. , o 2 , ... , 2 Vl Z l 



+-lim— -In 8 [3 + 3x l + 2xi+4:yiZi + -^-i-) . (3.9) 
2 i^oo 4' +i V X\ J 

The new variables xi, y\ and 2; fulfill recursion relations which are easy to deduce from 

their definitions and Eqs. 1)3.2113.5)1 . and are not difficult to iterate (starting with initial 

values X\ — j^, yi — 1 and Z\ — =, following from the numbers g% = 88, hi = 112, 

zi = 64 and ji = 112, found by direct computer enumeration of the corresponding HWs 

within the generator Gf(2)). One quickly finds that xi, yi, and Z\ tend to zero, and 

Zi <C xi <C yi, xf ~ for large Z, meaning that the last terms on the right-hand sides 

of equations ()3.8)) and (J3.9)) tend to zero. It is then straightforward to see that from 

(P and (jHll) h follows that 

,. Inflfj ln/ij In 
hm = hm = lim , 

and, according to 1)3.7)1 . one finds 

lncj = lim (3.10) 

Instead of the number j h which rapidly grows with /, it is convenient to introduce yet 
another variable 

Inji In 22/1 1\ 
' = 7; 5- 7 - 77 h (3- 11 ) 



4« 3 V4 A 1 

which has the initial value u\ = ln ^ 12 . Numerically iterating its recursion relation, 
together with those for xi, yi and zi, one can show that u\ tends to 1.2507788499... 
when I — > oo. Finally, since lncu = | lim/^oo u ; + ^In22, the connectivity constant for 
the 6 = 2 3d Sierpinski fractal is equal to u = 2.12587.... 
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In order to find the first correction to the leading-order behavior of C\ one needs to 
know the asymptotic behavior of the numbers xi, y\ and ji, according to the formula 
InCz+i In 16 \iaxi hiyi A ^ n 3i 

~A^+7 = iVm + iVm + iVm + iVm ' ( ] 

following from (J3.6|) and the definition of x\ and y\. Keeping only the leading-order terms 
in the recursion relation for x\ one gets xi + \ ~ constx 2 , which means that X\ behaves as 

xi ~ A 2 ' (3.13) 

for large I. Numerically iterating the recursion relations for xi, yi and z\ one finds 
A = lim^oo ^ = 0.9055.... It then follows that y l+1 w ^x x y h z l+1 fa ^Xiz h implying 
that & ~ 2 l , which, together with the numerically established relation xf ~ y^;, gives 

z x ~ 2- , / 2 x, ~ 2-'/ 2 A 2; , ^ ~ 2</ 2 A 2 ' . (3.14) 

On the other hand, from the definition of U\ (13. lip and the corresponding recursion 
relation it follows 

1 , / 4 „ 1 4 24 12 2 4 2 2 

= ui + }^ In 11 + — x fc + — x fc + —x k y k z k + — x fc y fc z fc + — y k z k 

In 22 /l 1 

+ 



3 V4 A 1 , 
Then, using (|3.1U|) . one can write 
In lln22 

;- = 21nW r 

4' 4' 3 

~ 1 , / 4 o 1 4 24 12 2 4 2 2 

~ 2^ ^ftt ln ( 1 + Yl Xk + Ti Xk + Yi XkVkZk + Yi Xk VkZk + \i Vk Zk 



and consequently, since 

~ 1 . / 4 3 1 4 24 12 2 4 

x ^fpt ln 1 1 + ji Xk + Yi Xk + YY Xfc2/feZfc + TY Xfc VkZk + ii 



fc=; 



, 104 " 1 11 104 

< ln Y —- = — t ln 

11 f-j4 fc+1 3 4' 11 



which is not difficult to show, one obtains 

hi ji = 2*4'lnu; + 0(l). (3.15) 
Finally, from (|3~T2J) - (f3~T5T) it follows that 

ln Q = Nilnu + JV, 1/2 In A^ + - ln N t + 0(1) , 

which means that the behavior (jl.l)) of the number of HWs, expected for homogeneous 
lattices, is also satisfied for this fractal lattice, with the following values of the exponents: 

1 1 
a = - , a = - . 

2 2 

The number of all possible HW configurations within a generator of the 3d 
Sieprinski fractal grows rapidly with b. In |Appendix C we give the corresponding 
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(b) 

Figure 4. Schematic representation of types of open Hamiltonian walks through (a) 
5-simplex and (b) 6-simplex of order I. 

recursion relations for the numbers gi, hi, i\ and ji found for the b = 3 3d SF, together 
with their initial values, and the formula for the number Ci + \ of closed HWs. The CPU 
time required for the enumeration and classification of HW configurations was so long 
for the b = 3 case that we could not go beyond it. However, the method used for b = 2 
3d SF in principle could be applied for any b with no qualitative difference. Analyzing 
recursion relations given in |Appendix C| in that manner, for the b = 3 3d SF we found 
that the number C\ of closed HWs obeys a scaling form similar to that of the b = 2 3d 
case with the following values for the connectivity constant uj and exponents a and a: 

In 3 

uj = 2.2722364... , a = , a = 0.386... . 

In 10 

4. Hamiltonian walks on n— simplex fractals 

To complete our analysis of the nature of Hamiltonian walks on 2d and 3d Sierpinski 
fractal lattices, we now turn to n-simplexes, which are in some ways a generalization 
of SFs for b = 2 in n — 1 dimensions. To obtain an n-simplex lattice one starts 
with a complete graph of n points and replaces each of these points by a new complete 
graph of n points. The subsequent stages are constructed self-similarly, by repeating 
this procedure. After I such iterations one obtains an n-simplex of order I, whereas the 
complete n-simplex lattice is obtained in the limit I — > oo. It is trivial to see that the 
connectivity constant uj for HWs on 3-simplex lattice is equal to 1, whereas Bradley 
found that uj = 1.399710... for the 4-simplex [THj. In principle, it is possible to establish 
an exact set of recursion relations for the numbers of suitably chosen prerequisite HWs 
on any ra-simplex, as Bradley did for n = 4. Here we'll demonstrate the method on 
n = 5 and n = 6 cases. 
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4-1. 5- simplex 

Any closed HW of order (/ + 1) on this fractal can be decomposed into 5 open HWs 
through 5-simplices of order I. There are two possible types of these open HWs, as 
depicted in Figure HJa). The first corresponds to walks which enter the simplex at one 
corner, visit all vertices inside it, and leave it - we shall denote the number of these walks 
by Cx t i- A walk of the second type enters the simplex at one of its five corners, wanders 
around it visiting some of the vertices inside it, leaves it through the second corner, 
and afterwards enters it again at the third corner, visits all the remaining vertices, and 
finally leaves it - let the total number of these walks on the Z-th order 5-simplex be C 2j ;. 
The total number Cj+i of closed HWs is equal to 

Q +1 = 12CI + 30CJA,, + 60C*! 3 A 2 / + 132C£, (4.1) 
which we found by computer enumeration, together with the recursion relations 

C 1>1+1 = 6C£, + 30C£,C2,i + 78(7^, + 96(7^ <7 2 3 ,* + 132<7 y <7 2 4 ,* + 132CJ, , 

C 2tl+1 = 2C\ X + 13(7^(7^ + 32Cj,Cf tl + 88(7^(7^ + 220<7 M <7 2 4 ,* + 186CJ, . (4.2) 

Initial values for these numbers are C\,\ = 6 and 6*2,1 = 2. From the recursion relation 
fl4~T|) it follows that 

= ±§* + JL (132 + 60*? + 304 + 12*f) , (4.3) 

where the new variable x\ = C\j,jCi \ satisfies the recursion relation 

_ 22 + 22xz + 16xf + 13xf + 5xf + xf 4 
Xl+1 ~ 186 + 220 xi + 88 xf + 32 xf + 13 xj + 2 xf ^ ' ^ 

obtained from (|4.2|) . Numerically iterating this relation, starting with x\ = 3, we find 
that xi — > 0.802318837... when I — > 00, and consequently, since Ni = 5 l , it follows that 

In a; = lim — ttt— = lim — -r- 2 - + lim — rrv In f 1 + — + —xf H xf^) = lim — . 

1^00 5 m 1^00 5 1 i->oc 5 m V 11 22 11 V 1^00 5 1 

The last limiting value can be quickly calculated if we introduce the variable 

lnC 2 ,i In 186 fl 1 \ In 2 

Vl ~ ~1T 4~ \5 ~ h l ) ~ T~ ' 

which, according to (|4.2|) . obeys the recursion relation 

1 , / 110 44 2 16 3 13 4 1 5 \ 

Vi+i = V/ + — ttt In H H xf H xf H xf H xf , (4.5) 

y + y 5 l+1 V 93 93 ' 93 1 186 1 93 V V ; 

and has the initial value y\ = 0. Then, iterating ()4.5|) simultaneously with ()4.4|) we find 

Vl -> 0.141065489481... for large/, and finally, since In a; = Iinii_ foo j/,+liil86/20+ln2/5, 

we obtain = 1.717769.... 

To examine the leading order correction to the number of Hamiltonian walks 

on the 5-simplex fractal, we note that for any k the quantity y\ can be written as 

Vl = Vk + YJm=k(Vm+i ~ Vm) ■ Taking the / — > 00 limit and keeping k fixed in that 

equation, one obtains: 



ln<7 2jfc In 186 ™ 1 / 110 44 2 16 3 13 4 1 5 
+ T^W + <^ k 5^W V "93" Xm + 93 Xm + 93 Xm + 186 Xm + 93^ 
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Since x\ = 3 and the array x m is monotonically decreasing, the sum on the right hand 
side of the above equation is bounded from above by 

00 1 111 1 

M2L72) E ^TT = ^(21.72)^(1 + - + ...) = - ^(21.72)^. 

m=k 

Therefore we can conclude that 

\nC 2 , k = 5 fc lncu + 0(l) and \nC hk = lnC 2jfe + \nx k = 5 k Into + 0(1). (4.6) 

Finally, for the number C\ of closed HWs, from ()4.3j) . we obtain the same behavior: 

In G = 5'lncj + 0(l). (4.7) 

Comparing with equation ()1.1|) we see that in the case of the 5-simplex fractal both the 
surface and the power correction to the number of HWs are absent. 



4-2. 6-simplex 

In addition to the Cij and C^z-type walks, already defined for the 5-simplex, one should 
introduce one other type of walks for complete enumeration of HWs on the 6-simplex 
lattice. These walks enter the 6-simplex three times, as depicted in Figure EJb). Let's 
denote their number on the 6-simplex of order / by C^. The total number Ci + \ of closed 
HWs within the 6-simplex of order (/ + 1) is equal to 

C l+1 = 60 Cj, + 360 C\f 2 .i + 1170 C\p\ x + 1920 C\f\ x + 3960 C^C^ 

+7920 CjjCZj + 5580 , (4.8) 

as we found by computer enumeration. Numbers Cij, C21 and C^i satisfy closed set of 
recursion relations, which can be numerically analyzed in a way similar, although more 
complicated than that used in the case of the 5-simplex lattice (see |Appendix D ). In 
contrast to the 5-simplex case, where numbers Cij, Ci,\ and C\ have the same asymptotic 
behaviour ()4.6|) and (|4.7p . here we obtained 

\nC hl = 6'lnw + 2 m lnA + 0(l), (4.9) 
\nC 2 ,i =Q l \noo + 2 l \n\ + 0(l) , (4.10) 
lnC 3ii =6 l \ncu + 0(l), (4.11) 

with u = 2.0550..., A = 0.9864 and 

In a =Nt\noj + N?lnns + 0(l), (4.12) 

with Ni = 6 l , a = In 2/ In 6 and fis — The surface correction to the number of HWs 
is present here, but still there is no term proportional to IniVi, which would imply the 
presence of the power correction. The same conclusion was obtained by Bradley jjlS.j for 
the 4-simplex lattice, with a — 1/2. 
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3 4 5 6 7 0.0 0.1 0.2 0.3 0.4 0.5 

Z 1/b 

Figure 5. (a) Connectivity constant li for Hamiltonian walks on Sierpinski fractals 
(y for 2d SF, and full triangles for 3d SF) and n-simplex lattices (0) as functions 
of the coordination number z together with previously found results for 4-simplex 
[13], hexagonal [12] (star), square [29] (■), triangular [30] (A), and cubic [15] (full 
diamond) lattices, as well as Flory [28] and Orland [21] approximations, top = (z — l)/e 
and uio — z/e, respectively, (b) Connectivity constant for 2d SFs as function of the 
reciprocal of the scaling parameter 6, where line connecting triangles serves merely as 
the guide to the eye. Open up-oriented triangle on vertical axes depicts the value of u 
for triangular lattice. 

5. Discussion and conclusions 

The values of the connectivity constant uj found by us, as well as some values previously 
found by other authors [T3] are depicted in Figure a), as functions of the coordination 
number z\ of the lattice, together with the Flory arid Orland |2I| approximations, 
and uj for hexagonal [T2|, square |30J, triangular [3T] and cubic [T3] lattices. One can 
clearly see that uj increases with z, which is in accord with Flory ujp = (z — l)/e and 
Orland et al too = z/e formulas, but it is obvious that uj depends on other lattice 
properties too. Furthermore, all values of uj lie between the values predicted by these 
two formulas, and it seems that these kind of approximations give satisfying results for 
fractal lattices studied here. Based on the information in this figure one may conclude 
that ujo and up are good upper and lower bounds for fractal HW connectivity constants. 
In Figure Efb) the connectivity constant uj for 2d SFs is presented as a function 

| The coordination number z (the average number of nearest neighbours per site) for the n-simplex 
lattice is equal to n, whereas it can be shown that z = 6(6 + 2)/(b + 4) for 2d SF with the scaling 
parameter b, z = 6 for b = 2 3d SF and z = 6.75 for b = 3 3d SF. 
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of the reciprocal of the fractal scaling parameter b. It seems that for b — > oo the 
connectivity constant might approach its triangular lattice value. This is not surprising, 
since for b = oo the first step of the construction of the corresponding 2d SF is already 
the wedge of the triangular lattice. In the same limit, 3d SFs approach the corresponding 
three-dimensional Euclidean lattice, so it would be useful to obtain recursion relations 
for HWs on these fractals for larger b, since there are fewer results for more realistic 
three-dimensional lattices. This task requires faster computers, as well as establishing a 
better algorithm for enumerating the HW configurations within a generator of a fractal, 
which is something we are planning to do in the nearest future. 

As for the correction to the leading-order asymptotic behavior of the number of 
HWs, we have shown that the surface correction fi^ appears for neither the two- 
dimensional Sierpinski fractals nor for the 5-simplex lattice. In contrast, for both three- 
dimensional Sierpinski fractals considered here, as well as for the 4-simplex ^3] and 
the 6-simplex lattice, the surface correction is present, with the value of the exponent 
a = l/df, where df is the fractal dimension of the corresponding lattice§. This result 
is certainly not a simple generalization of the formula proposed for the homogeneous 
lattices: o = (d — l)/d. This is actually not surprising: The correction term fig" in 
(II. 1J) was originally introduced in order to take into account possible surface tension 
effects, since at low temperatures a SAW forms a compact globule (see, for instance, 
|36|). The value a = (d — l)/d for homogeneous lattices then follows from the fact that 
the surface of such a globule is proportional to N^ d ~ 1 ^ d . In the case of fractal lattices 
it is, however, questionable whether such surface effects exist at all. For instance, all 
sites of the b = 2 3d SF lie on the surface, which is not the case for the b = 3 3d SF, 
but for both lattices the number of HWs has the correction term fj,^" . Or, in the case 
of an n-simplex lattice all sites have the same number n of neighbors, and one should 
not expect any surface correction, but still, for some of them the correction fif 7 was 
found. So, it seems that the term , obtained for some of the studied fractals, does 

not originate from the surface effects, which is also supported by the fact that N l l df is 
proportional to the mean radius of the globule formed by a HW, and not to its surface. 

It is interesting to note here that the existence of the term n$ in the scaling 
form for the number of HWs, coincides with the existence of the polymer coil to globule 
transition on the corresponding lattice. As was shown earlier, a polymer chain in a 
solvent, modeled by SAWs on the 2d Sierpinski fractals [27] and 5-simplex |32J can 
exist only in the swollen phase; on the 3d Sierpinski fractals jSHl EH ES] , 4-simplex |2S] 
and 6-simplex [32] lattices, when the temperature is lowered the polymer undergoes a 
collapse transition from an expanded state to a globule state (compact or semi-compact 
[3*4*j). Analyzing the asymptotic behaviour of different types of open HWs on fractals 
one can observe that the collapse transition exists on lattices whose topology allows for 
the statistical domination of HWs which are not localised. In particular, HWs on the 
2d SG fractals cannot enter a generator of any order more than once, i.e. all the walks 

§ Fractal dimension for 3d SF is df — ln(6(6 + l)(b + 2)/6)/ln6, and for the n-simplex lattice 
df = Inn/ In 2. 
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are localised. On the 5-simplex fractal this is possible (see Figure Ufa)), but the number 
C\ } i of localised HWs and the number of delocalised HWs (walks that enter every 
/th order 5-simplex twice) are of the same order (see ()4.6j) . i.e. the delocalised HWs do 
not dominate. On the other hand, the delocalised HWs on the 3d SG fractals (j-type, 
see Figure 0} and on the 6-simplex (C^jr and C^-type, see Figure Hfb)), as well as on 
the 4-simplex are possible, and furthermore, the number of these walks is much 
larger than the number of the localised HWs (see section |3] and subsection I4.2j) . This 
observation strongly resembles the conclusion obtained in a series of recent papers [21] 
about the delocalisation of knots in the low-temperature globular phase. Although the 
term 'delocalisation' was not used in quite the same sense in these two contexts, it seems 
that the same effect is in question, and this problem deserves further investigation. 

Finally, the power dependence of the overall number Cat of closed HWs on the 
number of sites N of the lattice was found only for the 3d SFs. A more detailed 
inspection of the calculation of the exponent a for these lattices reveals that 

lnfi 

a = const lim — . 

l^oo I 

The numbers hi and i\ correspond to the localised open Hamiltonian configurations that 
visit the maximal (4) and the minimal (2) number of vertices, respectively, within the 
generator of order I (see FigureEJ). In the case of 2d SFs, numbers of open configurations 
visiting the maximal or the minimal number of vertices were hi and gi, respectively, and it 
was shown that the ratio hi/gi=const for every I. Consequently, lim^oo \n{hi/ii)/l = 0, 
which may be the formal explanation for the absence of the power correction to the 
number of HWs on 2d SF. On the other hand, on the n-simplex lattices only one type of 
localised HW configurations is possible, so it appears that the power term in the scaling 
form for the number of HWs is obtained on lattices where a larger number of different 
types of localised configurations is possible. 

In conclusion we can say that the method of exact recursion relations turned out to 
be very powerful for the generation and the enumeration of extremely long Hamiltonian 
walks on the two- and three-dimensional Sierpinski and n-simplex fractals. Furthermore, 
it allows for a detailed numerical analysis of HWs of different topologies. This enabled 
us to find various scaling forms for the number of closed HWs on these lattices. In the 
case of two-dimensional Sierpinski fractals, a closed-form expression is obtained for the 
connectivity constant. Very interesting results were obtained for the three-dimensional 
Sierpinski fractals. This should be utilized for attaining deeper insight into the realistic 
physical problems which can be modeled by Hamiltonian walks. 
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Figure Al. h— and g— type HW on a Gf, x generator, which can be obtained one 
from another in a unique way. These two walks differ only in the parts surrounded by 
circles. 



No: 1634) and JS from the United States Department of Energy under contract to the 
University of California. 

Appendix A. Recursion relations for HWs on 2d Sierpinski fractals 

In this appendix we present the derivation of relations ()2.1j) and ()2.3j) for open and 
closed HWs on 2d SFs, for general b. 

It is obvious that the number of closed HWs on the (/ + l)th order generator 
is of the form Ci + \ = J^iBih^gf", where is the number of all closed Hamiltonian 
configurations consisting of «j steps of /i-type and $ steps of g-type ("step" is here a part 
of HW that traverses /th order generator within the considered (/ + l)th order generator). 
By definition, closed HW visits all sites of the (/ + l)th generator, consequently it 
traverses all b(b + l)/2 /th order generators within it, so that 

+ = (A.l) 

On the other hand, every h-type step occupies all three vertices of the /th generator it 
traverses, whereas a g-type step occupies only two vertices (the entering and the exiting 
one). This means that the numbers ctj and (3i also have to satisfy the following equation 

2a , +ft= (^±M±^ (A . 2) 

since the number of /th order vertices inside the (/ + l)th generator is (6+l)(6+2)/2. The 
system of equations (jA.l|) and (jA.2|) has a unique solution «j = 6+1, /3, = (b+l)(b— 2)/2 
which completes the derivation of the relation ()2.1j) . 

For open /i-type HWs, in a similar way, we have hi + \ = J^iAih^g^, where 
exponents Xi and yi satisfy the system Xi + i/i = b ( b +^ , 2xi +yi— ( fc+1 ^ 6+2 ) — 1 ; whose 
only solution is Xi = b, = b(b — l)/2. At the same time, every g-type HW on a 
generator Gf +1 can be obtained from one and only one h-type HW (and vice versa), by 
substituting one h-step with a g-step, as depicted in Figure Al. This means that all of 
these walks have exactly (b — 1) /i-steps and [b(b — l)/2 + 1] g-steps, i.e. relations (j2.3J) 
are correct. 
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Appendix B. RG equations for SAW model on 2d SF 

In this appendix we give exact RG equations (|2.1()|) for calculating the connectivity 
constant for the ordinary SAW model on the 6 = 4 and the 6 = 5 two-dimensional 
Sierpinski fractals. These equations were obtained via computer enumeration and 
classification of all SAW configurations within the corresponding fractal generator. RG 
relations (j2.10|) for 6 = 6 and 7 were also found, but they are too cumbersome to be 
quoted here and they are available upon request. 



6 = 4: 

B' = B 4 + QB 5 + QB 6 + 9B 7 + 9B 8 + 9B 9 + AB 10 + AB 3 B 1 + 30B 4 B 1 + 36B 5 B 1 
+56B 6 B 1 + 58B 7 B 1 + 56B 8 B 1 + 2QB 9 B 1 + QB 2 B 2 + 57 B 3 B 2 + 8AB 4 B\ 
+ 13AB 5 B 2 + U3B 6 Bl + 128S 7 S? + 56B 8 B 2 + ABB 3 + 51B 2 B 3 + 96B 3 B 3 
+ 15QB 4 B 3 + 168B 5 B 3 + 128B 6 B 3 + A0B 7 B 3 + B\ + 2\BB\ + 55B 2 B\ 
+93B 3 B 4 + 9AB 4 B 4 + A8B 5 B 4 + 3B\ + \ABB\ + 28B 2 B\ 
+ 205 3 J B 1 5 + B\ + 455? 
5; = B 6 B 1 + 65 7 5i + 7B S B 1 + 45 9 5i + 65 5 5 2 + 38B 6 B 2 + 445 7 5 2 + 265 8 5 2 
+1AB 4 B 3 + 925 5 5 3 + 1025 6 5 3 + 5QB 7 B 3 + 165 3 5 4 + 10QB 4 Bf + 10AB 5 B 4 
+ A0B 6 B 4 + 9B 2 B\ + 58B 3 B\ + 405 4 5 1 5 + 255? + \2B 2 B\ 
6 = 5: 

B' = B 5 + 105 6 + 205 7 + 305 8 + 545 9 + 685 10 + 985 11 + 945 12 + 86B 13 
+38B 14 + 165 15 + 5B 4 B 1 + 60B 5 B 1 + U0B 6 B 1 + 228B 7 B 1 + U3B 8 B 1 
+586B 9 B 1 + 867B 10 B 1 + 8545 n 5i + 786B 12 B L + 3485 13 5i + U0B 14 B 1 
+ 10B 3 B 2 + U6B 4 B 2 + A02B 5 B 2 + 718B 6 B 2 + \52\B 7 B\ + 2Y37B 8 B\ 
+3203B 9 B 2 + 32405 10 5 2 + 29185 n 5 2 + 12685 12 5f + 4585 13 5 2 + 10B 2 B 3 
+ 18AB 3 B 3 + 6105 4 5 3 + 1218S 5 S? + 28A6B e B 3 + A316B 7 B 3 + 64335 8 5 3 
+66485 9 5 3 + 56305 10 5 3 + 23065 n 5 3 + 6645 12 5 3 + 5BB 4 + 126B 2 B 4 
+523B 3 B 4 + 12095 4 5 4 + 31705 5 5 4 + 53075 6 5 4 + 76785 7 5 4 + 79605 8 5 4 
+59605 9 5 1 4 + 21045 10 5 4 + 3605 14 5 4 + 5? + AABB\ + 2A9B 2 B\ + 7105 3 5 1 5 
+21595 4 5 1 5 + 41185 5 5 1 5 + 5604S 6 J B 1 5 + 55545 7 5 1 5 + 32925 8 5 1 5 + 7765 9 5 1 5 
+65? + 5955? + 2345 2 5? + SgiE 3 ^ 6 + 20315 4 J B 1 6 + 2A79B 5 Bf + 2086 J B 6 J B 1 6 
+744fi 7 J B 1 6 + 5B 7 + 365.B 7 + 21452.B 7 + 630 J B 3 J B 1 7 + 62QB 4 B 7 + 32AB 5 B 7 
+ B 8 + 28BB 8 + lir^ 2 ^ 8 + 72 J B 3 J B 1 8 + 2B\ + 

5; = B 8 B 1 + 12_B 9 5! + 39B 10 B 1 + 48fi 11 J B 1 + GOfi 12 ^ + 3AB 13 B 1 + lGB 14 ^ 
+8B 7 B\ + 102 J B 8 J B 1 2 + 3AAB 9 B 2 + A32B W B 2 + 556B U B 2 + 3\AB 12 B\ 
+ 1A0B 13 B 2 + 27B 6 B 3 + Seefi 7 ^ 3 + 1278B 8 B 3 + 1616B 9 B 3 + 2098B W B 3 
+ 115QB U B 3 + A58B 12 B 3 + 50B 5 B 4 + 722B 6 B 4 + 2600B 7 B 4 + 32545 8 J B 1 4 
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+4128fi 9 J B 1 4 + 2128B 10 B* + 66AB u Bf + ttB 4 B\ + 852B 5 B 5 X + 3U8B & B\ 
+3808B 7 B* + AA7AB S B\ + l968B 9 Bf + 360ff 10 £ 1 5 + 36B 3 B® + 6lOB 4 Bf 
+23025 5 J B 1 6 + 2590£ 6 £ 1 + 25A0B 7 Bf + 7365 8 J B 1 6 + YiB 2 B\ + 2545 3 J B 1 7 
+981 J B 4 J B[ + 9A8B 5 B 7 + 592B 6 B 7 + 255 s + 545 2 J B 8 + 2205 3 J B 8 

+ \AAB 4 B\ + ABBf + 20B 2 B? . (B.l) 

Appendix C. Recursion relations for HWs within the b = 3 3d SF 

Here we give the recursion relations for the numbers gi, hi, %\ and ji of open HW 
configurations within the b = 3 3d SF, obtained via computer enumeration. With 
symbols g l+1 , h i+1 , g h h h ji and i x abbreviated to g', ti, i' , f, g, h, i and j, 

respectively, these relations have the following form: 

g' = 6120i 2 g 3 j 5 + 3312i 2 g 4 j 4 + 8176i 2 g 5 j 3 + 5068i 2 g 6 j 2 + 296Ai 2 g 7 j + 776i 2 g 8 
+13296i 3 gfh + 12688i 3 g 2 j 4 h + 36832i 3 g 3 fh + 3650Ai 3 g 4 j 2 h + 27768i 3 g 5 jh 
+9200i 3 g 6 h + 2080i 4 j 4 h 2 + 20224i 4 #j 3 /i 2 + A0A6Ai 4 g 2 j 2 h 2 + A606Ai 4 g 3 jh 2 
+ 21856i 4 g 4 h 2 + 28A8i 5 j 2 h 3 + 12192i 5 gjh 3 + 11328i 5 g 2 h 3 + 512i 6 h 4 

ti = 2928ig 4 f + I288ig 5 j 4 + 3296ig 6 j 3 + 1760ig 7 j 2 + 936ig s j + 232ig 9 + I3296i 2 g 2 fh 
+10608i 2 g 3 j 4 h + 28160i 2 g 4 j 3 h + 238A0i 2 g 5 j 2 h + 16632i 2 g 6 jh + 5168i 2 g 7 h 
+12768i 3 fh 2 + U336i 3 gj 4 h 2 +50176i 3 g 2 j 3 h 2 + 63232i 3 g 3 j 2 h 2 + 5662Ai 3 g 4 jh 2 
+21856i 3 g 5 h 2 + 1382Ai 4 j 3 h 3 + 3AA32i 4 gj 2 h 3 + A8928i 4 g 2 jh 3 + 27968i 4 g 3 h 3 
+ 582Ai 5 jh 4 + 8192i 5 gh 4 

i! = 1250Ai 3 g 2 f + 7880i 3 g 3 j 4 + 192A0i 3 g 4 j 3 + 13232i 3 g b j 2 + 822Ai 3 g 6 j + 218Ai 3 g 7 

+528i 4 j 5 h + A592i 4 g] 4 h+2A22Ai 4 g 2 fh + 36928iVj 2 /i + 33728i 4 g 4 ]h+\329Qi 4 g b h 
+118Ai b fh 2 + 838Ai 5 gj 2 h 2 + 18880i 5 ^ 2 j/i 2 +13024i 5 ^ 3 /i 2 + QA0i 6 jh 3 + lQ00i 6 gh 3 

f = A308i 2 f + 5808i 2 gj 7 + 17A2Ai 2 g 2 f + ll93Qi 2 g 3 f + 19lQAi 2 g 4 j 4 +U09Qi 2 g 5 f 
+9208i 2 g 6 j 2 + 3360iVj + 544i 2 # 8 + llQlQi 3 fh+2lAA0i 3 gfh + 5102Ai 3 g 2 j 4 h 
+66096iVj 3 /i + 56056i 3 ^ 4 j 2 /i+28864i 3 ^ 5 j/i + 6400i V>i + 10312i 4 j 4 /i 2 
+35296i 4 ^ 3 /i 2 +53248i 4 ^ 2 j 2 /i 2 + A5AA0i 4 g 3 jh 2 + U080i 4 g 4 h 2 + 5728i 5 j 2 h 3 
+ l25AAr>gjh 3 + 71Q8i 5 g 2 h 3 + 512i 6 h 4 . 

The initial values of these numbers are g\ = 497000, h\ = 728480, i\ = 340476, 
ji = 811468 and the formula for the number Cj+i of closed HWs within Gf +1 (3) is: 

C l+1 = 92gfjf + ASgfjt + 8g}° + 1792t l9 f h l3 f + V2A8i l g 7 h l2 i + 38Ai l gfh l 

+75Q8t 2 g 4 h 2 j 2 + 710At 2 gfh 2 ji + W08i 2 g?h 2 + 10560^^ 2 /i 3 j' 2 + 134402 3 ^ 3 / i 3 j / 
+7680i 3 gfh 3 + AA80i 4 h 4 jf + 7Q8Qi 4 g l h 4 l3l + 6016i 4 gfh 4 + 512^/if . 
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Appendix D. Analysis of recursion relations for HWs on 6-simplex lattice 

Recursion relations for the numbers Cu, C 2 , and C 3j ; of open HW configurations within 
the 6-simplex of order I have the following form: 

C[ = 5544C 2 C 4 + i728C 3 Cf C 3 + 2592C 3 C 3 + 120C?Cg + 480C*C 2 C 3 + 960C?Cf 

+48C 1 5 C 3 + 216C 1 5 C 2 + 24C? + 25008C 4 C 3 2 + 20544C 2 5 C 3 + 6576C 2 6 + 11328C 1 C 2 3 C 3 2 

+15264CiC , 2 4 C 3 + 8688CiC 2 5 +4992C 2 C 3 C 3 (D.l) 

C 2 = 94336C 2 C 4 + 76800C 3 C 3 + 48160C 4 C 3 2 + 23520C 2 5 C 3 + 6576CiC 2 5 + 17120CiC 4 C 3 

+16672CiC 3 C 3 2 + 2832C 1 2 C 2 C 3 2 + 5088C 2 C 3 C 3 + 832C 3 C 2 C 3 + 3620C , 1 2 C 4 

+1232C 3 C 3 +144C 4 C 2 C 3 + 324C 1 4 C 2 +64C 1 5 C 2 + dCf + 16C?C 3 (D.2) 

C 3 = 541568C 3 6 + 94336C 2 3 C| + 43200C 4 C 3 2 + 14448C 2 5 C 3 + 2940C 2 6 + 6252C 1 C 4 C 3 

+2568C 1 C 2 5 + 1416C 2 C 3 C 3 + 954C 1 2 C 4 + 208C 3 C 2 3 +54C 4 C 2 2 

+6C 1 5 C 3 + Y1C\C 2 + Cf , (D.3) 

where G[ = Cy+i and Cj = C^/, with the initial values C^i = 24, C 2 .\ = 6, C 3j i = 1. 
Introducing new variables 

Cu C 2 i lnC 3/ In 541568 (\ 1 



X ' C 2 / y ' C 3 / 24 6' 5 V6 6' 

one can obtain closed set of recursion relations which iterates towards xi, y\ — > 0, z\ 
0.280204... for I — ► 00. On the other hand, from fl£g|) it follows that 
lnQ + i _ lnyj ^ In 5580 1 / 44 22 2 32 3 13 4 2 5 1 6 
6 l+i ~ 6 i + 6 *+i + 6 w ln v i + 31 X/+ 31 Xi + 93 X ' + 62 X/ + 31 X; + 93 X ' 
whereas from ()D.2|) one gets 

lnC 2 , ;+1 _ lnC 3 ,i ln yi In 94336 , 1 / 600 1505 2 521 2 
6 m " qi + V 1 + 6 i+i + 6 i+i [n { 1 + 737^ + 2948^ + 2948 + 
1^ „2„,2 1 ^35 3 535 3 159 2 3 1 5 4 3 6 4 



-xtyf H u, H xiyf H xfy, H x, y, H x,y, 

5896 ' ' 2948 1 2948 w 2948 ' 1 1474 47168 ,w 

and, consequently, 

lnC 3) / 4 lnj/f 
in to = hm : 1 — hm — y- . 

1^00 6 3 1^00 6 

Since the numerical analysis shows that lnyi/6 l 0, we finally obtain u = 2.0550... 

To find the leading order asymptotic behavior of the number of Hamiltonian walks 

on the 6-simplex fractal (Cy,C2j, C 3j / and Ci), we conduct an analysis similar to the 

one used in the case of the 5-simplex. However, we can tell right away that Ci t i,C 2> i and 

c c 

C31 will have mutually different asymptotics, since their ratios x\ = ^ and yi = ^ 
were found to be approaching zero for large I. By iterating the recursion relations for 
x\ and yi, one finds that the ratio of X[ and yi quickly approaches a constant equal to 
1.521868... Keeping only the terms with the lowest sum of powers in x\ and y\ in the 
recursion relations, we find Xi + i = 12 , yi + \ = 2^^^-, for I ^> 1 and the ratio 
%1+i/yi+i is indeed 1.521868... From the above equation for yi+i, one can see that 

In yi m 2 l In A, A = const = lim ^ = 0.9864... (D.4) 

l— >oo 2 
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The asymptotic relation (|4.11|) can be obtained starting with z\ = Zk+J2m=k( z m+i — 
z m ) and following a procedure completely analogous to the one used in the 5-simplex 
case. Relations ()4.9j) and (|4.1()j) then follow from ()4.11|) . (|D.4|) and the fact that xi and 
yi are proportional in the large I limit. From (|4.9|) - ([4.11|) it is apparent that in the 
limiting case C\ t i <C Cij, <C C^,u since A < 1. Therefore, it holds that Ci + \ ~ 5580C| ; , 
so finally In Ci = 6 l \nu + 3 * 2 l In A, which is equivalent to (|4.12j) . 
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